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Abstract 

Angle Resolved Photoemission Spectroscopy (ARPES) measurements on high temperature 
superconductors, such as BiSrCaCuO, show three main components. These are a quasiparticle 
spectral peak that develops below the superconducting transition temperature Tc, an 
accompanying broad background of secondary electrons and a dip feature beside the main 
quasiparticle peak. The broad background may originate from inelastic processes in which the 
photoelectron emits and absorbs spin fluctuations. Calculations of the quasiparticle spectral 
weight are presented incorporating these spin fluctuation mediated inelastic processes in which 
the development of the superconducting gap A has been incorporated into the magnetic 
susceptibihty xiq^E). A dip feature develops below Tc in the quasiparticle spectral weight due 
to the shifting of spin fluctuation spectral weight, lnix{q,E), from low energies to energies 
greater than 2A. These results predict that the dip feature in the ARPES spectrum in high 
temperature superconductors such as BiSrCaCuO, is evidence for the opening of a spin gap 
below Tr. 



Introduction 

Angle Resolved Photoemission Spectroscopy (ARPES) has developed into an important 
probe of the superconducting and normal state properties of the cuprate-oxide, high temperature 
superconductors. ARPES measurements support LDA predictions for the electronic 
bandstructure of the cuprates [|l| and a (1^2 _y2 symmetry superconducting order parameter [|[. 
ARPES provides information on how this gap evolves with doping and its connection to the 
normal state gap seen in underdoped cuprates [^. 

The overall shape of the ARPES spectra measured in these experiments is not fully 
understood, however. A common feature of the data is the co-existence of a broad background 
and a peak. The latter is interpeted as the quasiparticle spectral weight peak of the Fermi surface 
electrons emitted from the crystal after the absorption of the ultraviolet photon (~ 20ey). The 
relative heights of the background and the peak vary from one experiment to another [0, H, 0]. In 
addition to the background, a dip feature is seen to develop when the crystal is cooled oelow the 
superconducting transition temperature Q. The connection between this dip feature and a similar 
feature observed in tunneling experiments on the cuprates [^, 0, § has been the subject of recent 
investigations by the present authors |^, |lO| . The tunneling density of states is another measure of 
the quasiparticle spectral weight. The difference between tunneling and ARPES arises, in part, 
because tunneling measurements provide an average of the quasiparticle spectral weight along a 



line of states in the Brillouin Zone, determined by the directional tunneling matrix element |1C]. 
ARPES measurements provide information about the quasiparticle spectral weight in a small 

region, determine d by experimental resolution, around a fc-space point on the Fermi surface. 

Recent work [10| proposed a common explanation for the broad ARPES background and the 



linearly increasing tunneling conductances seen in many high temperature superconductor tunnel 



junctions |11]. It was proposed that the ARPES background is due to the simultaneous emission 
or absorption of spin fluctuations by the electron as it is escapes from the surface layer of the 
crystal after absorbing the ultraviolet photon. The same processes can occur in the surface layer 
near the tunnel barrier in a tunnel junction and can lead to a wide variety of tunneling 



conductances. Related work has been carried out by others 11^ , 13 1. 

Previous work [^] made use of a phenomenological model for the spin fluctuation spectral 
weight which had been used to flt inelastic neutron scattering measurements of the spin 



fluctuation spectrum in YBaCuO |14]. This model for the spin fluctuation spectrum did not have 



superconducting correlations incorporated in it. This issue has become relevant due to the 
observation of structure in the spin fluctuation spectrum which appears to be directly connected 



to the development of superconductivity in the cuprates |15|. Another important issue is the 
relative height of the main spectral peak and the background as well as its dependence on the 
underlying magnitude of the spin fluctuation susceptibility. The overall magnitude of the inelastic 
background contribution to the quasiparticle spectral weight was multiplied by a constant fltting 
factor in pO|] in order to reproduce results that compare favorably with ARPES data. 

In the present work, the spin fluctuation spectral weight, D{E), which determines the 
inelastic background, is calculated using the RPA approximation for the underlying electronic 
spin susceptibility, xiQ^^)- The present calculations incorporate the tight binding bandstructure 
with next nearest neighbor hopping, appropriate for materials of interest such as BSCCO, and 
superconductivity arising from a d^2__y2 symmetry order parameter. The effects of the 
superconducting state on the spin fluctuation spectrum are incorporated into a calculation of 
A{k, E) and a comparison is made between predictions and experimental ARPES measurements 

of A{k^ E) above and below Tc- 

One of the results of this work is the prediction that the dip feature seen in ARPES data 
below Tc on BSCCO B is indirect evidence of the development of the spin-gap in x(^, E) due to 
the onset of supercondfuctivity. The present work allows a controlled investigation to be carried 
out of the dependence of the relative magnitude of the inelastic background and the main 
quasiparticle peak on the strength of the electron-spin fluctuations. Finally, the Van Hove peak in 
the underlying tight binding bandstructure which can strongly influence the energy dependence of 



the spin fluctuation spectral weight D(E), is included. 

The present model for the quasiparticle spectral weight A{k, E) assumes that the main peak 
in the ARPES spectrum arises from photoelectrons created by the absorption of ultraviolet 
photons by electrons at the Fermi energy. The accompanying broad background is due to the 
simultaneous absorption or emission of spin fluctuations by other electrons as they absorb 
photons. Spin fluctuation emission dominates at the low temperatures of interest here. In the 

ARPES experimental technique, the momentum or k space region being probed can be identifled 
with relatively good precision using energy and momentum conservation. The photoelectrons 
which contribute to the inelastic background are secondary electrons, which having emitted spin 

fluctuations, are scattered into the same momentum direction, k, as those electrons which yield 
the main ARPES spectral peak. They are collected by the detector and labelled with the same 

momentum k as the main peak yielding an overall APRES spectrum for a particular k vector. 

Theoretical Model 

The total quasiparticle spectral weight, A{k^ E), in our model is made up of two 
contributions. The main peak is calculated from 

A°{k,E) = --ImG{k,E) (1) 

TT 

where 

The electronic bandstructure is defined as $,k = —2t{cos{kx) + cos{ky)) — At cos{kx)cos{ky) — n 
[Bwhere we have chosen a typical value for the cuprates of t = — 0.45t. The superconducting 
order parameter is given by A^ = Ao{cos{kx) — cos{ky)/2 with Aq = O.lt. An electronic damping 
parameter is also included through T = To + TiiT /Tc)'^ with Fq and Fi chosen to be 0.04t and 
0.05t respectively. The chemical potential ^ is chosen to be — 1.75t. 
The inelastic contribution is calculated from 

A,nei{k,E) = / dEImG{k,E) (3) 

[D{E - E'){n{E - E) + f{-E'))<S){E - E) 
+ D{e' -E){n{E' -E) + f{E'))@{E' -E)] 

where n{E) denotes the Bose-Einstein distribution and f{E) denotes the Fermi function. 

The spin fluctuation spectral weight, D{E), is defined as the momentum integrated spin 
fluctuation susceptibility multiplied by the square of the electron-spin fluctuation coupling 
constant, g/t, 

D{E) = -i-2 |_ J dqx J_^ dqyig/tflmxiq, E) (4) 

where 



The imaginary part of the bare electronic spin susceptibihty is defined as 

\v^x\q. E) = Sp. / dE {f{E + E) - f{E )) (6) 



.2 „.2 



[{[Upj^qUp + UpVpUp+qVpj^q)]ImG{p + Q, E + E)ImG{p, E ) 
+ [ul+qvl - UpVpUp+qVp+q]ImG{p + q, E + E)ImG{p, -E ) 
+ [^l+q^l - UpVpUp+qVp+q]ImG{p + q,-E - E)ImG{p, E ) 
+ [Vp+q'Vp + UpVpUp+qVp+q]ImG{p + q,-E - E)ImG{p, -E )] 

where Up and Vp represent the usual superconducting coherence factors and 

In.ai^,E) = l ^^J^^,^^, (7) 

where Ep = \/{ip + A^). The real part of x'^io'^E) is obtained by Kramers-Kronig transform. 
The total quasiparticle spectral weight is then given by 

A{k, E) = A"{k, E) + aA,nei{k, E) (8) 

where a is an overall constant multiplicative factor which determines the relative contributions of 

the elastic A°[k, E) channel and the inelastic channel Aineiik, E). For the results presented in this 
paper, a is chosen to be in the range from 1 to 4, depending on the values chosen for the other 

physical parameters in the calculation of A{k,E). The choice of a incorporates the relative 
weighting of the bulk (leading to A°(k,E)) and surface physics (leading to Ainei{k,E)). 

In our previous [1^] work based on |1^, the expression for the total spectral weight of 
equations (3) and (8) was extracted from the expression for the total elastic and inelastic 
tunneling current which can be written as 

I{V) = j dENs{E)NN{E - eV)[f{E) - f{E - eV)] (9) 



Ns{E) = SAT^'f{A"{k,E)+aA,nei{k,E)) (10) 



where 

NsiE) = L^ 

In obtaining equations (3) and (8), the tunneling matrix element squared in the inelastic 
contribution to the current, I(V), (denoted by |A(^'-^)p in equation (19.37b) of ||l^) was assumed 
to be given by jA^^'^^p = a\A^^'^>\'^(g/t)'^. By replacing ei with e/, — eV in equations (19.33b) and 
(19.37c) of [^], the combined elastic and inelastic current can then be written as in equation (9) 
(using the notation |T|'p in place of lA'-^-'^'p). The spectral weight function of equation (8) yields 
the density of states Ns{E) using equation (10). 

The inelastic tunneling channel described here implies that an electron tunnels into the 
superconducting crystal along a direction determined by the directional tunneling matrix element 
T^'' and then emits a spin fluctuation which brings in the spin fluctuation coupling constant g/t 
into the overall tunneling matrix element. The spectral weight function of equation (8), which 
determines the density of states in equation (10), can then be used to interpet ARPES data. 

One approximation in the present approach is the lack of conservation of momentum. This 



approximation is also a feature of related work in this field |12, |l^]. This approximation is valid 
for the case of sufficiently high electronic disorder and scattering which will broaden the 

underyling A°{k,E). In the clean limit, equation (3) should contain terms involving 



T,g'ImG{k + q,E )Imx{q, E — E ) within the integral on the right hand side of the equation which 
would require a significant increase in numerical computation to yield an accurate answer for 
Ainei{k,E). The assumption inherent in the present work is that ImG{k + q,E ) is sufficiently 
broadened by disorder that it is a reasonable approximation to take it outside the sum over spin 
fluctuation wave vectors g and replace it with ImG{k,E ). This approximation can also be 
further justified by noting that Imxiq, E — E ) is strongly peaked at Q = (vr, vr) and that for the 

electronic wave vectors of interest here ImG{k, E ) ~ ImG{k + Q,E ). We have also tested that 
the results for the total quasiparticle spectral weight to be presented in the accompanying figures 
can be generated with higher values of damping F, and correspondingly more broadened A"{k,E), 
than have been used in the work shown here. 

Results 

Figures (1) and (2) show results for the spin fluctuation spectral weight D{E) for the case of 
a constant coupling constant g = U. Figure (1) depicts D{E)t at T /Tc = 1.0 and 0.3 The effect of 
the onset of superconductivity is evident in the figure in the removal of spectral weight from 
below 2Ao to higher energies. 

In the figures depicting D{E)t in the present work, the underlying value of 'Eq'Imxiq, E) is 

given by dividing D{E)t by (U/t)'^. For example, the peak value for D{E)t of approximately 0.3 

in figure (1) implies an underlying value of 'Sqimxiq, E) of 2.0 states per eV assuming a value of 

t = ISOmey. 

In calculating the two curves in figure (1), the sum rule 

1 /^a /•+°°^^ Imx{q,E) 

— 7T a q (IE- — — — — -— - = constant (11) 

27r2i ^i_oo {l-exp{-E/kBT)) ^ ' 

is imposed. The same value of [/ = l.Ot was used for the two D{E) curves in figure (1) with 
/i = -1.75t. 

Figure (2) depicts D{E)t for the case U = 1.5t which is close to the largest possible value of 
U/t for the choice of ;U = —1.75t before the RPA approximation for xiqi £')breaks down. 

The resulting quasiparticle spectral weight curves, A{k,E)t, for these D{E)t curves are 

depicted in figures (3), (4) and (5). All A{k,E)t curves presented in this paper are calculated for 

k on the Fermi surface at k = (vr, 0.1624). Figure (3) depicts A{k,E)t for the D{E)t of figure (1), 
with Q = 4.0 in equation (8). The development of the dip feature associated with the onset of the 
spin-gap in D{E)t is clearly visible below the superconducting transition temperature. 

Figures (4) and (5) depict A{k,E)t using the D{E)t of figure (2) and use a = 1 for figure (4) 
and a = 1.5 for figure (5). The large increase in magnitude of D{E)t in figure (2) compared to 
that of figure (1) allows a sizeable inelastic background to occur from equation (8) with the elastic 
and inelastic channels contributing about equally. Unlike [|l^], the connection between a and the 
magnitude of the underlying Imx{q,E), which is determined by U/t from equation (5) can now 
be explored in a controlled manner. The effect of interactions in D{E)t is to shift the peak 
slightly below 2Ao as can be seen by comparing D{E) in figures (1) and (2). This effect can result 

in a disappearance of the dip feature in A{k,E)t as can be seen in figure (4). However, a slight 
increase in a restores the feature as shown in figure (5). Figure (6) is generated using the D{E)t 
of figure (2) with a = 3.0 and Fq = 0.06t and Fi = 0.075t in equation (2). These figures illustrate 
the wide variation in the relative heights of the main elastic peak and inelastic background in 
A(k, E)t that can be generated within the present model. A similar variation is seen in ARPES 
experiments suggesting that the background is not a bulk phenomenon but instead a reflection of 
the surface physics which can probably vary from sample to sample. 



The peak at E = 2Ao in D{E)t in figures (1) and (2) is caused by a strong peak in the 
underlying Imx{q,E) at q= (7r,7r) at ^ = 2Ao for ^ = — 1.75t. This peak shifts downwards in E 
as U increases due to the part of the denominator involving RexiQ, E) in equation (5) for the 
susceptibility. As has been just pointed out, the extent to which this occurs can influence the 
ability to produce a dip feature in the spectral weight. The peak &t E = 2Ao is also sensitive to 
the choice of chemical potential //. A slightly smaller negative value for // will result in the peak 
moving to higher energies in Imx{q,E) and diminishing in height |jl^. The E = 2Ao peak in 
Imx{q, E) is the result of the underlying Van Hove peak in the tight binding bandstructure and 
moves to different values of E as the chemical potential is varied. 

Figures (7) and (8) depict D{E)t for the choice g = —Jo{cos{kx) + cos{ky)) |[l^ which 
enhances the role of the (vr,-7r) peak in Imxiq,E). Jo= l.Ot and 1.2t for figure (7) and Jo = 1.2t 
and 1.4t for figure (8). This choice of values for Jq ensures that the sum rule of equation (11) is 
satisfied H^. The resulting A{k,E)t are depicted in figures (9) (which uses the D{E)t of figure 
(7)) and (10) (which uses the D{E)t of figure (8)) for a = 4.0 and a = 2.0 respectively. In 
calculating Aineiik,E)t from equation (3), the spectral weight D{E)t is integrated over energy E 
and, as a result, sharp peak structure in D[E)t is somewhat smeared out in the resulting A{k^ E)t. 

The effect on A{k,E)t of choosing a different point in k space has been investigated before in 

pO[| . Anisotropy as a function of k enters into the calculation of the quasiparticle spectral weight 
in several ways. Apart from the underlying electronic bandstructure S,k, the the most significant 
sources of anisotropy are in the order parameter A^ and in strong coupling effects such as the 
quasiparticle damping rate, T. This last issue was treated phenomenologically in |jl^ by 

increasing the value of T for those regions of k space where A^ = 0. The overall effect is to smear 
out the main quasiparticle peak and, as a result, eliminate the dip feature in the quasiparticle 
spectral weight. 

The spin fluctuation spectral weight D{E)t can also be used to estimate the quasiparticle 
damping rate, denoted by T in equation (2), from 

r = 2vr / dE . /;^7; (12) 

JO sinn[E/ kB-L ) 

Using the D{E)t from figure 1(b) for T/Tc = 1 and assuming Aq = 2kBTc, T is found to be 0.06t. 
This value is comparable to values of the quasiparticle damping rate F used in generating the 
results of figures (1) to (10). 

Equation (12) is derived from the conventional strong coupling treatment of spin 
fluctuations. The model for the contribution of the inelastic background in this paper is different 
to conventional treatments of strong coupling effects in calculations of the quasiparticle spectral 
weight and tunneling densities of states. The role of spin fluctuations in the cuprates has been 
widely investigated in both the normal and superconducting states|2C, 21, ^]. The resulting 



quasiparticle spectral weights A{k, E) in these calculations can be used to generate density of 

states curves from N{E) = T,k\Tk\'^A{k,E). In the superconducting state, densities of states 
curves will display small corrections relative to the underlying weak coupling densities of states. 
This will not provide an explanation for the rapidly increasing tunneling densities of states at high 



bias voltages measured in tunneling experiments on the cuprates [11]. This type of variation with 
bias voltage is a signature of an additional inelastic channel, which in the cuprates is assumed to 
involve the emission and absorption of spin fluctuations in the surface region of the sample. 
The difference between the inelastic tunneling model and conventional strong coupling 
approaches can also be seen by considering the type of Feynman diagrams that arise in the usual 
calculation of the tunneling current using linear response theory. Conventional strong coupling 
corrections are incorporated with diagrams of the type shown in figure (103) combined with figure 
(100a) of [|l6|. Inelastic tunneling is calculated using a different diagram as depicted in figure (98) 
of [jl6[ where both vertices of the diagram are joined by the propagator representing the spin 
fluctuation. 



Conclusions 

The spin gap below E = 2Ao in Imx{q-,E) in figures (1), (2), (7) and (8) occurs because a 
spin fluctuation must have an energy E greater than or equal to this threshold in order to create 
a quasiparticle-quasihole pair at low temperatures. The effects of this on the quasiparticle 
spectral weight and tunneling density of states have been investigated before in different ways 



1 23, 24 1 . [23 1 incorporated this type of pairbreaking physics phenomenologically in a model based 
on the marginal Fermi liquid theory for s-wave superconductivity and investigated the resulting 
quasiparticle spectral weight and density of states curves. In [|j], a quasiparticle damping 
mechanism based on the same pairbreaking mechanism for a 6.^2 _y2 order parameter was 
incorporated approximately into the quasiparticle spectral weight and the resulting S-I-S current 
characteristics were calculated for comparison with experiment Q. 

The work of [|^, which is the basis for the present calculations, used a phenomenological 
model for the spin fluctuation spectral weight. In [^], a dip feature is present in the quasiparticle 
spectral weight due to a combination of the narrowing of the main peak, because of the reduction 
of the scattering rate in the superconducting state, and the underlying shape of the model used for 
Imx{q,E). No superconducting correlations were incorporated into the model for xiQ^E) which 
are now known to be important [^] and the dip feature was an accidental feature of the model. 

In conclusion, results for a model of the quasiparticle spectral weight A{k, E)t have been 
presented which incorporate a conventional elastic peak and an inelastic background arising from 
spin fluctuation emission processes; equation(3). The goal is to interpet ARPES measurements on 
high temperature superconductors. The results of the present work provide a model connecting a 
microscopic calculation of the spin fluctuation susceptibility xiQ,E), equations (5) and (6), to the 
magnitude and overall shape of the inelastic background seen in ARPES. The spectral weight 
curves generated in this approach can also be used to interpet tunneling conductance 
measurements on high temperature superconductors [|lO[. The dip feature seen in some ARPES 
data is caused by the development of a spin gap in the underlying spin susceptibility at the onset 
of superconductivity. 
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Figure Captions 
All horizontal axes in the figures are in units of E/t. 

Figure (1): The monientum integrated spin fluctuation spectral weight D{E)t from 
equation (4) for g = U = l.Ot. T/Tc = 1.0 (solid line). T/Tc = 0.3 (line with symbols). 

Figure (2): The momentum integrated spin fluctuation spectral weight D{E)t from 
equation (4) iov g = U = IM for T/T^ = 0.3. 

Figure (3): The quasiparticle spectral weight A{k,E)t for k on the Fermi surface, a of 
equation (8) is 4.0. This curve is generated using D{E)t of figure (1). k is chosen to be (vr, 0.1624) 
on the Fermi surface. T/Tc = 1-0 (solid line). T/Tc = 0.3 (line with symbols). Curve 3(b) depicts 
the contributions of the elastic and inelastic channels separately for T/Tc = 0.3. The contributions 
of the two channels depicted in figure 3(b) are typical of the other A{k, E)t curves in this paper. 

Figure (4): The quasiparticle spectral weight A(k,E)t for k on the Fermi surface, a of 

equation (8) is 1.0. This curve is generated using D{E)t of figure (2). k is chosen to be (vr, 0.1624) 
on the Fermi surface. 

Figure (5): The quasiparticle spectral weight A(k,E)t for k on the Fermi surface, a of 

equation (8) is 1.5. This curve is generated using D{E)t of figure (2). k is chosen to be (vr, 0.1624) 
on the Fermi surface. 

Figure 6: The quasiparticle spectral weight A{k, E)t for k on the Fermi surface, a of 
equation (8) is 3.0. This curve is generated using D{E)t of figure (2). Fq = 0.06t and Fi = 0.075t 
in equation (2). k is chosen to be (tt, 0.1624) on the Fermi surface. 

Figure (7): The momentum integrated spin fluctuation spectral weight D{E)t from 
equation (4) for g = Jq = —Jo{cos{qx + cos{qy)). T/Tc = 1-0, Jo = 1.2t (dashed line). T/Tc = 0.3, 
Jo = l.Ot (line with symbols). 

Figure (8): The momentum integrated spin fluctuation spectral weight D(E)t from 
equation (4) for g = Jq = —Jo{cos{qx + cos{qy)). T/Tc = 1-0, Jo = 1.4t (solid line). T/Tc = 0.3, 
Jo = 1.2t (line with symbols). 

Figure (9): The quasiparticle spectral weight A{k,E)t for k on the Fermi surface, a of 

equation (8) is 4.0. This curve is generated using D{E)t of figure (7). k is chosen to be (tt, 0.1624) 
on the Fermi surface. T/Tc = 1-0 (solid line). T/Tc = 0.3 (line with symbols). 

Figure (10): The quasiparticle spectral weight A{k,E)t for k on the Fermi surface, a of 
equation (8) is 2.0. This curve is generated using D{E)t of figure (8). k is chosen to be (tt, 0.1624) 
on the Fermi surface. T/Tc = 1-0 (dashed line). T/Tc = 0.3 (line with symbols). 
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